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1.0  Scope  and  Objectives 


Giannotti  and  Buck  Associates,  Inc.  has  been  tasked  by  the 
Naval  Ship  Engineering  Center,  Code  6136,  to  perform  motion  analy- 
ses on  an  oceanographic  and  research  survey  ship,  T-AGOS,  now  in 
the  concept  design  stage.  Of  particular  interest  in  the  study  are 
two  points  on  the  ship  representing  equipment  interface  locations. 
The  first  point  is  located  at  Frame  59  on  ship  with  three  alternate 
heights  having  been  investigated.  The  second  point  is  located  on 
the  upper  side  of  the  cap  rail  on  the  fantail  bulwark.  The  pri- 
mary vehicle  used  in  the  analysis  is  the  Five-Degree-of -Freedom 
Ship  Motions  Computer  Program  maintained  by  NAVSEC  Code  6136. 
Statistical  characteristics  of  the  motions  and  long  term  motion 
predictions  were  desired.  This  report  presents  the  approach  to 
the  problem  and  reports  the  results  of  the  analyses. 

\ 

2.0  The  Five-Degree-of-Freedom  Ship  Motions  Program 

The  basic  5-DOF  Program  predicts  ship  response  in  five  degrees 
of  freedom  (surge  neglected)  and  was  developed  under  the  super- 
vision of  Dr.  C.  Chryssostomidis  of  MIT  using  the  theory  presented 
by  Salvesen,  Tuck  and  Faltinsen  (1).  For  the  purpose  of  calculat- 
ing the  various  hydrodynamic  coefficients  used  in  the  program  the 
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ship  sections  ere  represented  by  so  called  Lewis-forms.  These 
are  ship-like  sectional  shapes,  completely  fair,  without  chines, 
knuckles  or  discontinuities  and  are  obtained  by  a conformal 
mapping  technique  using  local  section  beam,  draft  and  area  as 
arguments.  The  difficulties  Introduced  in  representing  a double 
chine  hull  with  this  form  will  be  discussed  in  the  next  section. 

The  program  contains  several  options  for  seaway  representa- 
tion. In  addition  to  regular  wave  responses  at  any  heading  in 
waves  of  any  length,  heights,  or  slope,  long  crested  unidirec- 
tional irregular  waves  may  be  introduced  having  spectral  charac- 
teristics either  of  the  standard  Pierson-Moskowitz  spectrum  or  of 
any  desired  input  spectrum.  The  Pierson-Moskowitz  spectrum  is 
usually  used  for  representation  of  a fully  developed  storm  seas  of 
long  duration  but  such  seas  occur  less  frequently  in  nature  than 
do  partially  developed  seas  of  shorter  duration.  Seas  of  this  type 
are  better  represented  by  the  Bretsnelder  two  parameter  spectrum. 

As  a third  alternative, wave  spectra  based  on  direct  observation 
rather  than  mathematical  modeling  can  be  introduced. 

For  the  purposes  of  this  study  the  Pierson-Moskowitz  spectrum 
was  used  throughout. 

The  program  computes  absolute  and  relative  vertical  motions 
amplitudes,  velocities,  and  accelerations  and  transverse  accelera- 
tions only  at  any  point  specified - 

In  its  basic  form  the  program  also  predicts  hull  bending  moments, 
torsional  moments  and  shear  forces  generated  dynamically  in  waves. 
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In  addition,  the  occurrence  and  frequency  of  deck  wetness,  slamming, 
and  propeller  racing  may  be  forecast  although  none  of  these  fea- 
tures were  used  in  this  study. 

A modification  to  the  original  program  permits  the  introduc- 
tion of  experimentally  determined  Response  Amplitude  Operators,  a 
feature  which  was  Invoked  in  this  study. 

2.1  Program  Modifications 

An  experimental  study  of  several  round  bilge  and  hard  chine 
hull  forms,  including  the  T-AGOS  parent  hull  form  has  been  con- 
ducted by  Slbul  (2)  at  the  University  of  California.  Sibul  com- 
pared the  resistances  of  the  hulls  in  calm  water  and  their  motions 
in  head  and  beam  seas.  Past  experience  has  Indicated  that  analy- 
tical ship  motions  programs  give  generally  good  results  in  pre- 
dicting the  smaller  motions-pitch,  heave,  yaw  and  sway — but  rolling 
motion,  because  of  its  usually  large  amplitude  and  non-linear 
character,  is  more  difficult.  In  this  case  the  difficulty  is 
compounded  by  the  fact  that  the  hard  chines  are  not  well  modeled 
by  the  Lewis-form  representation.  As  Sibul  showed  experimentally 
the  hard  chines  have  a substantial  effect  in  damping  rolling 
motions. 

The  5-DOF  Program  contains  a feature  which  permits  the 
simultatlon  of  bilge  keels  in  the  representation.  Since  these 
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appendages  have  a generally  similar  effect  to  the  sharp  corners 
of  the  hard  chine  hull  some  initial  computer  runs  were  made  in  an 
attempt  to  duplicate  the  experimental  evidence  using  bilge  keels 
to  simulate  the  chines.  Figure  1 shows  the  results  of  several  of 
these  attempts  compared  to  the  Sibul  data.  The  ordinate  of  these 
curves  is  roll  amplitude  divided  by  wave  slope.  As  the  waves 
become  very  long,  that  is  as  the  frequency  goes  to  zero,  the 
response  should  be  Just  equal  to  the  forcing  function.  In  other 
words,  the  non-dlmensicnalized  rolling  motion  should  approach 
unity.  As  indicated  in  Figure  1 the  modified  5-DOF  Program  does  not 
predict  this,  a fact  which  casts  a shadow  on  the  value  of  the  roll 
predictions  made  by  the  program.  The  fault  could  be  either  in  the 
theory  or  in  the  way  the  theory  is  programmed.  It  will  take 
further  study  to  determine  the  precise  cause  and  correct  it.  In 
the  meantime,  efforts  to  use  the  analytical  roll  prediction  for 
this  application  were  abandoned  and  the  Sibul  experimental  roll 
operators  were  introduced  in  the  program.  Analytically  calculated 
operators  in  the  other  modes  of  motion  were  continued. 

In  Table  1 Sibul 's  experimental  results  are  reproduced  together 
with  the  predictions  made  by  the  modified  5-DOF  Program.  As  Table  I 
shows,  the  modified  Program  using  analytically  determined  operators  in 
pitch  and  heave  and  experimentally  determined  operators  in  roll 
slightly  overpredicts  in  pitch,  underpredicts  in  heave  and  under- 
predicts in  roll.  Although  the  agreement  is  not  impressive  it  is 
reasonable  considering  differences  in  the  experimental  and  theo- 
retical input  wave  spectrum. 
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TABLE  I 


Comparison  Between  Experimental  Results  (Ref.  2)  and 
the  Modified  5-DOF  Motion^  Program  Prediction 


Mode 

Heading 

Spee-i,  kts. 

* — SIBUL-*  ^-MODIFIED  5 
Sig.  Dbl.  Sig. 

Wave  ht.  Amp.  Wave  Ht. 

DOF  — * 
Dbl. 

Amp. 

Pitch,  deg. 

Head  Seas 

3.0 

8.64 

7.64 

8.64 

7.87(1) 

Heave , ft. 

Head  Seas 

3.0 

8.64 

6.19 

8.64 

4.44(1) 

Roll,  deg. 

Beam  Seas 

0.0 

9.60 

31.8 

10.20 

26.5(2) 

(1)  Using  analytically  predicted  RAO  with  1.0'  bilge  keel  width. 

(2)  Using  RAO  determined  experimentally  in  regular  wave  tests. 


2 . 2 Angular  Motions 


Motions  at  specified  points  are  of  special  interest  in  this 
study  in  that  these  points  are  candidate  equipment  interface 
locations.  The  coordinate  system  and  the  locations  of  the 
specified  points  are  shown  below. 


tijo  9*0  Wo  Vo  Wo  s*o  Vo  j'o  z'o  Vo  a 


The  origin  of  the  (x,y,z)  coordinate  system  is  fixed  with  respect 
to  the  mean  position  of  the  ship  in  the  undisturbed  free  surface 
with  the  z-axis  vertically  upward  through  the  center  of  gravity. 

The  translational  displacements  in  the  x,y,z  directions  are  r^. 
rt2»  and  r^.  The  rotational  displacements  about  the  x,y,  and  z axes 
are  n^,  Hj.  and  n6» 
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Thus 


- surge 

n4  " 

roll 

n2  * sway 

n5  " 

pitch 

r)^  ■ heave 

n6  " 

yaw 

The  locations  of  the  specified  points 

are: 

of  interest 

X 

y 

z 

Point  1 -5.84 

0.0 

64.00 

Point  2 -5.84 

0.0 

70.63 

Point  3 -5.84 

0.0 

77.25 

Point  4 -102.84 

0.0 

9.25 

Points  1,  2,  and  3 are  three  alternate  heights  of  a sensor  to  be 
located  at  Frame  59.  Point  4 is  located  on  the  cap  rail  of  the 
fantall  bulwark. 

As  noted  the  program  compute  s absolute  and  relative  vertical 
motions,  amplitudes,  velocities  and  accelerations  and  transverse 
accelerations  only  at  specified  points.  In  the  basic  MIT  program 
this  feature  computed  correctly;  however,  in  the  version  that 
accommodates  externally  generated  RAO's  there  is  an  error  in  the 
program  which  has  not  as  of  this  time  been  identified  although 
the  other  features  of  the  program  compute  correctly.  Considering 
the  nature  of  the  equipment  involved,  the  vertical  and  transverse 
motions  are  primarily  of  Interest  at  Point  4 and  for  this  purpose 
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the  predictions  generated  by  the  basic  program  are  considered 
satisfactory. 

For  the  equipment  to  be  located  at  Points  1,  2,  or  3 the 
angular  displacements,  maximum  angular  velocities,  maximum  angular 
accelerations,  and  maximum  x,y,z  components  of  linear  velocity  and 
acceleration  at  the  location  of  the  point  are  of  greater  interest. 

An  attempt  was  made  to  Introduce  computation  of  these  elements  Into 
the  external  RAO  version  of  the  program.  Within  the  time  and  cost 
constraints  of  this  study  the  attempt  was  not  successful  and  was 
deferred  to  a later  effort.  To  produce  useful  and  timely  answers 
hand  calculation  was  resorted  to  for  a limited  number  of  the  most 
significant  cases.  Since  rolling  motions  are  of  largest  magnitude 
the  maximum  values  in  the  roll  mode  of  angular  velocity,  angular 
acceleration,  tangential  velocity  (linear  y-component  at  the  point), 
tangential  acceleration  (linear  y-component  at  the  point),  and 
normal  acceleration  (linear  z-component  at  the  point)  were  com- 
puted by  hand.  These  values  were  also  computed  for  the  worst  case 
in  the  pitch  mode  to  give  a comparison  of  magnitudes. 

The  computational  relationships  involved  in  calculating  the 
angular  motions  are  given  below. 

Maximum  angular  velocity: 

n.  ■ u>  n.  . i ■ 4,5,6 

1 el 

max  max 


Maximum  angular  acceleration: 


max 


Maximum  tangential  velocity:  v£  - 

i 

Maximum  tangential  acceleration:  at  ■ R^f^ 
Maximum  normal  acceleration: 


max 


Here, 


*"  ■ Rx  rtW)2 

' R*'r'imax'1^™ax 

toe  ■ frequency  of  encounter 

- M(1^  co.  a 

s 

vg  ■ ship  speed,  ft/sec 
a ■ heading  angle 

Rk  " projected  radius  from  ship  c.g.  to  ktl1  point 
in  plane  normal  to  motion  axis . 


In  our  case  Rk  - + OG  - + 5.58' 

Point  1^  - 69.58 
Point  2:R2  - 76.21 
Point  3:R3  - 82.83 

Note  that  if  values  of  tangential  velocity  and  acceleration  and 
normal  acceleration  are  given  for  one  point,  values  at  other 
points  may  be  determined  by  simple  ratio. 


3.0  Results  of  Short  Term  Predictions 


The  predictions  cited  in  this  section  are  all  for  irregular 
long  crested  seas  using  the  Pierson-Moskowitz  Spectrum.  Significant 
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wave  heights  used  to  represent  the  various  Sea  States  are  shown 
In  Figure  2 which  Is  a graphical  display  of  the  traditional  table 
of  Wind  and  Sea  Scales  due  to  Wilbur  Marks. 

Three  headings — head  seas  (180°) , bow  seas  (135°)  and  beam 
seas  (90°),  and  three  speeds,  0,  3,  and  5 knots — were  addressed 
In  the  study. 

Roll  Is  the  dominant  mode  of  motion  In  the  seakeeping 
performance  of  the  T-ACOS  design,  and  the  mode  which  primarily 
drives  the  motions  of  Points  1,  2,  and  3.  Figure  3 illustrates 
the  interaction  of  the  ship  in  the  roll  mode  and  the  fully  devel- 
oped seaways.  The  usual  operator  approach  is  employed: 

Seaway  Energy  Spectrum,  S^(u>)  x Ship  Response  Amplitude  Operator, 

Y (w)  ■ Ship  Response  Spectrum,  S_  (u>) . 

th  n** 

Note  the  sharply  tuned  character  of  the  Roll  RAO.  Because  of 
this  feature  the  ship  does  not  respond  to  much  of  the  huge  blocks 
of  energy  In  the  State  6 and  State  7 Sea  States.  Any  reduction 
in  the  Metacentrlc  Height  (GM)  during  the  design  process  will 
cause  the  ship  to  roll  at  a lower  natural  frequency  thus  shift- 
ing the  RAO  to  the  left  thus  producing  much  larger  rolling  motions 
and  much  larger  motions  at  Points  1,  2,  and  3. 

The  significant  heaving,  pitching,  and  rolling  motions,  that 
is  the  average  of  the  1/3  highest  motions,  are  displayed  in 
Figures  4,  5,  and  6.  Heaving  motions  will  drive  the  motion  of 
Point  4 on  the  fantail.  Figure  4 shows  that  these  motions  are 
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strongly  dependent  on  Significant  Wave  Height,  but  only  weakly 
dependent  on  speed  and  heading.  Weak  dependency  on  speed  is 
also  exhibited  in  the  pitching  and  rolling  motions  displayed 
in  Figures  5 and  6.  The  fact  that  the  roll  response  is  flaten- 
lng  out  with  increasing  wave  height  is  due  to  the  relationship 
of  the  ship  natural  roll  period  to  the  period  of  maximum  energy 
in  the  Seaway  Spectrum  as  discussed  above  and  depicted  in  Figure  3. 

Figures  7,  8,  and  9 tell  the  story  of  the  motions  of  Points 
1,  2,  and  3 in  the  roll  mode  which  as  been  noted  will  be  the 
dominant  mode.  Values  for  Point  3 only  are  displayed  but  values 
for  points  1 and  2 may  be  inferred  by  applying  the  ratios  indi- 
cated on  the  figures.  As  would  be  expected  the  trends  are 
similar  to  those  shown  in  Figure  6.  Velocities  and  accelerations 
of  these  points  don't  really  get  much  worse  beyond  a State  6, 
but  then  State  6 is  a severe  seaway  for  a 204-ft.  ship.  Figure  9 
compared  to  Figure  8 shows  that  normal  accelerations  are  much 
smaller  than  the  tangential  accelerations. 

Figure  10  presents  the  predictions  for  vertical  motions  at 
the  fantail  point,  Point  4.  The  amplitudes,  velocities  and 
accelerations  are  all  quite  moderate  and  again,  show  a rather 
weak  dependency  on  heading  and  speed. 

The  foregoing  data  are  tabulated  in  Table  II  which  gives 
the  ship  motion  predictions  and  Table  III  which  gives  the  point 
motion  predictions.  Table  III  also  gives  one  case  of  point 
motion  prediction  in  the  pitch  mode  from  which  it  may  be  Inferred 
that  point  motions  from  this  source  are  modest  compared  to  those 
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Note:  To  obtain  tangential 
velocities  at  Points  1 & 2 
multiply  Point  X values  by 
.840  and  .920  respectively. 
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SIGNIFICANT  HEIGHT 


Fig.  8 
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Fig.  10  TRENDS  OF  VERTI  CAL  MOTIONS 


H D G.  135° 


“ HDG.  180° 

C/5 


ACCELERATION 


t-J  5 ns. 

^ } 3 K T S. 

2^- — zi “}  o ns. 


POINT  4 


HDG.  90°.  0,  3,  5 RTS, 


10  12 


22  24  26  28  30  32 


SIGNIFICANT  HEIGHT,  FT. 


stemming  from  roll.  Transverse  accelerations  at  Point  4 are 


also  Included. 

It  should  be  borne  in  mind  that  all  of  the  above  deal  with 
short  term  predictions,  that  is  motions  resulting  from  a storm 
condition  which  may  remain  stationary  for  a matter  of  only  a few 


hours.  In  the  next  section  we  will  build  on  this  data  base  to 


TABLE  II 


RESULTS  OF  SHORT  TERM  SHIP  MOTION  PREDICTIONS1 
(Average  highest  motions  are  given  In  double  amplitudes) 


STATE  4 

STATE  5 

STATE  6 

2 

Slg.Ht.  - 6.5  ft. 

Sig.Ht.  - 10.2  fg. 

Slg.Ht.  - 16.9  ft. 

Heading  Angle 

180°  135°  90 

180  135  90U 

180°  135°  90 

SPEED  0.0  kts . 

HEAVE 


RMS 

.59 

.86 

1.64 

1.45 

1.81 

2.60 

3.20 

3.58 

4.29 

Avg.1/3  highest 

2.38 

3.45 

6.56 

5.80 

7.24 

10.43 

12.80 

14.33 

17.16 

Avg.  1/10  highest 

3.03 

4.40 

8.37 

7.40 

9.23 

13.30 

16.31 

18.27 

21.88 

Avg.  1/1000  highest 

4.57 

6.64 

12.63 

11.17 

13.94 

20.08 

24.63 

27.58 

33.03 

PITCH 

RMS 

1.2 

1.3 

.2 

1.9 

1.8 

2 

2.8 

2.5 

. L 

Avg.1/3  highest 

4.7 

5.0 

.7 

7.7 

7.4 

.7 

11.3 

9.9 

.6 

Avg. 1/10  highest 

5.9 

6.4 

.9 

9.8 

9.4 

.9 

14.4 

12.6 

.8 

Avg.  1/1000  highest 

9.0 

9.7 

1.3 

14.8 

14.2 

1.3 

21.8 

19.0 

1.2 

SWAY 

RMS 

0. 

.47 

.94 

0. 

1.00 

1.6  9 

0. 

2.11 

3.25 

Avg.1/3  highest 

0. 

1.87 

3.77 

0. 

3.99 

6.78 

0. 

8.44 

12.99 

Avg. 1/10  highest 

0. 

2.39 

4.80 

0. 

5.08 

8.64 

0. 

10.76 

16.56 

Avg. 1/1000  highest 

0. 

3.60 

7.25 

0. 

7.68 

13.05 

0. 

16.24 

25.0 

ROLL 

RMS 

0. 

1.9 

3.0 

0. 

4.5 

6.6 

0. 

6.9 

10.0 

Avg.1/3  highest 

0. 

7.8 

12.1 

0. 

18.1 

26.5 

0. 

27.6 

40.1 

Avg. 1/10  highest 

0. 

9.9 

15.4 

0. 

23.0 

33.8 

0. 

35.2 

51.1 

Avg. 1/1000  highest 

0. 

15.0 

23.3 

0. 

34.8 

51.0 

0. 

53.2 

77.1 

YAW 

RMS 

0. 

.5 

.1 

0. 

.8 

.2 

0. 

1.2 

.3 

Avg.1/3  highest 

0. 

2.0 

.5 

0. 

3.2 

.8 

0. 

4.6 

1.1 

Avg. 1/10  highest 

0. 

2.6 

.7 

0. 

4.1 

1.1 

0. 

5.9 

1.4 

Avg. 1/1000  highest 

0. 

3.9 

1.0 

0. 

6.2 

1.6 

0. 

8.9 

2.1 

(Continued) 
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TABLE  II 
(Continued) 


RESULTS  OF  SHORT  TERM  SHIP  MOTION  PREDICTIONS1 
(Average  highest  motions  are  given  In  double  amplitudes) 


STATE  4 

STATE  5 

STATE  6 

2 

Heading  Angle 

Sig.Ht 

180° 

. - 6. 
135° 

5 ft. 

90° 

Sig.Ht.  - 10.2  ft. 
180°  135°  90° 

Sig.Ht.  ■ 16 
180°  135° 

.9  ft. 
90° 

Speed  3.0  kts. 

HEAVE 

RMS 

.65 

.95 

1.64 

1.54 

1.92 

2.61 

3.31 

3.70 

4.29 

Avg.1/3  highest 

2.59 

3.79 

6.56 

6.16 

7.67 

10.43 

13.24 

14.78 

17.15 

Avg.1/10  highest 

3.30 

4.83 

8.36 

7.86 

9.77 

13.29 

16.88 

18.85 

21.86 

Avg. 1/1000  highest 

5.0 

7.29 

12.62 

11.87 

14.76 

20,07 

25.48 

28.45 

33.01 

PITCH 

RMS 

1.2 

1.3 

.1 

2.1 

1.9 

.1 

3.0 

2.6 

.1 

Avg.1/3  highest 

4.9 

5.2 

.6 

8.2 

7.6 

.6 

12.0 

10.3 

.5 

Avg.1/10  highest 

6.2 

6.6 

.7 

10.5 

9.7 

.7 

15.3 

13.1 

.6 

Avg. 1/1000  highest 

9.4 

9.9 

1.1 

15.8 

14.7 

1.1 

23.1 

19.8 

.9 

SWAY 

RMS 

0. 

.43 

.94 

0. 

.92 

1.70 

0. 

1.95 

3.26 

Avg.1/3  highest 

0. 

1.74 

3.77 

0. 

3.68 

6.81 

0. 

7.80 

13.02 

Avg.1/10  highest 

0. 

2.21 

4.81 

0. 

4.70 

8.69 

0. 

9.94 

16.61 

Avg. 1/1000  highest 

0. 

3.34 

7.26 

0. 

7.09 

13.12 

0. 

15.0 

25.07 

ROLL 

RMS 

0. 

1.8 

3.0 

0. 

4.4 

6.6 

0. 

6.9 

10.0 

Avg.1/3  highest 

0. 

7.1 

12.1 

0. 

17.5 

26.5 

0. 

27.7 

40.1 

Avg.1/10  highest 

0. 

9.0 

15.4 

0. 

22.3 

33.8 

0. 

35.3 

51.1 

Avg. 1/1000  highest 

0. 

13.7 

23.3 

0. 

33.7 

51.0 

0. 

53.3 

77.1 

YAW 

RMS 

0. 

.4 

.1 

0. 

.7 

.2 

0. 

1.0 

.3 

Avg.1/3  highest 

0. 

1.8 

.5 

0. 

2.8 

.8 

0. 

4.1 

1.1 

Avg.1/10  highest 

0. 

2.2 

.7 

0. 

3.6 

1.0 

0. 

5.2 

1.4 

Avg. 1/1000  highest 

0. 

3.4 

1.0 

0. 

5.5 

1.6 

0. 

7.9 

2.0 
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TABLE  II 
(Continued) 


RESULTS  OF  SHORT  TERM  SHIP  MOTION  PREDICTIONS1 
(Average  highest  motions  are  given  In  double  amplitudes) 


STATE  4 STATE  5 STATE  6 


2 

Heading  Angle 

Slg.Ht. 

180° 

■ 6.5 
135° 

ft. 

90° 

.2  ft. 
90° 

Slg.Ht 

180° 

. ■ 16 . 
135° 

9 ft. 
90° 

Speed  5.0 

kts . 

HEAVE 

RMS 

.72 

1.03 

1.63 

1.64 

2.02 

2.60 

3.42 

3.80 

4.28 

Avg.1/3  highest 

2.87 

4.13 

6.52 

6.57 

8.08 

10.39 

13.68 

15.20 

17.13 

Avg.1/10  highest 

3.66 

5.26 

8.31 

8.38 

10.30 

13.25 

17.44 

19.38 

21.84 

Avg. 1/1000  highest 

5.52 

7.95 

12.55 

12.65 

15.56 

20.0 

26.34 

29.26 

32.98 

PITCH 

RMS 

1.2 

1.3 

.1 

2.1 

1.9 

.1 

3.1 

2.6 

.1 

Avg.1/3  highest 

4.9 

5.1 

.5 

8.4 

7.7 

.5 

12.3 

10.4 

.4 

Avg.1/10  highest 

6.3 

6.5 

.6 

10.8 

9.8 

.6 

15.7 

13.3 

.5 

Avg. 1/1000  highest 

9.5 

9.9 

1.0 

16.3 

14.9 

1.0 

23.8 

20.1 

.8 

SWAY 

RMS 

0. 

.41 

.95 

0. 

.88 

1.71 

0. 

1.85 

3.26 

Avg.1/3  highest 

0. 

1.65 

3.78 

0. 

3.51 

6.84 

0. 

7.41 

13.04 

Avg.1/10  highest 

0. 

2.10 

4.82 

0. 

4.47 

8.72 

0. 

9.45 

16.63 

Avg. 1/1000  highest 

0. 

3.18 

7.28 

0. 

6.75 

13.16 

0. 

14.27 

25.10 

ROLL 

RMS 

0. 

1.5 

3.0 

0. 

3.6 

6.6 

0. 

5.8 

10.0 

Avg.1/3  highest 

0. 

5.9 

12.1 

0. 

14.3 

26.5 

0. 

23.0 

40.1 

Avg. 1/10  highest 

0. 

7.5 

15.4 

0. 

18.2 

33.8 

0. 

29.4 

51.1 

Avg. 1/1000  highest 

0. 

11.3 

23.3 

0. 

27.5 

51.0 

0. 

44.3 

77.1 

YAW 

RMS 

0. 

.4 

.1 

0. 

.7 

.2 

0. 

1.0 

.3 

Avg.1/3  highest 

0. 

1.6 

.6 

0. 

2.7 

.9 

0. 

3.8 

1.3 

Avg. 1/10  highest 

0. 

2.1 

.7 

0. 

3.4 

1.2 

0. 

4.9 

1.6 

Avg. 1/1000  highest 

0. 

3.1 

1.1 

0. 

5.1 

1.8 

0. 

7.3 

2.5 

^Uslng  experimentally  generated  roll  RAO's 
Angle  between  direction  of  ship  and  direction  of  waves. 
180°  • head  seas,  135°  - bow  seas,  90°  ■ beam  seas. 
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TABLE  IH 


RESULTS  OF  SHORT  TERM  POINT  MOTION  PREDICTIONS 
(RMS  of  maximum  values) 


1 


STATE  4 

STATE  5 

STATE  6 

2 

Heading  Angle 

ft. 

90° 

Sig.Ht.  - 10.2  ft. 
180°  135°  90° 

Sig.Ht.  - 16.9  ft. 
180°  135°  90° 

POINT  33,  SPEED  0 kts. 


ROLL 


Ang.Vel. rad/sec  ^ 

0. 

0.0253 

0.0450  0. 

Ang. Acccl . rad /sec 

0. 

0.0194 

0.0505  0. 

Tang. Vel. ft/sec 

0. 

2.10 

3.72  0. 

Tang .Accel .ft/s ec2 

0. 

1.61 

4.18  0. 

Norm. Accel . ft /sec 2 

0. 

0.124 

0.27  0. 

POINT  3 

, SPEED  3 kts. 

ROLL 

Ang.Vel. rad/sec 

0. 

0.0365 

0.0450  0. 

Ang. Accel. rad /sec 2 

0. 

0.0250 

0.0505  0. 

Tang. Vel. ft/sec 

0. 

3.02 

3.72  0. 

Tang. Accel .ft/sec2 

0. 

2.07 

4.18  0. 

Norm. Accel. ft /sec2 

0. 

0.721 

0.27  0. 

POINT  3 

. SPEED  5 kts. 

PITCH 

Ang. Vel. rad/sec 

Ang .Accel . rad/sec2 
Tang. Vel. ft/sec2 
Tang. Accel .ft/sec2 

Norm. Accel. ft /sec 2 

0.0560 

0.0745 

0. 

0.0857 

0.1259 

0.0415 

0.0655 

0. 

0.0615 

0.0935 

4.69 

6.17 

0. 

7.10 

10.42 

3.44 

5.43 

0. 

5.09 

7.75 

0.663 

0.90 

0. 

1.51 

3.09 

0.0600 

0.0745 

0. 

0.0932 

. 

i 

0.1259 

0.0480 

0.0655 

0. 

0.0730 

0.0935 

4.97 

6.17 

0. 

7.72 

10.42 

3.98 

5.43 

0. 

6.05 

7.75 

0.770 

0.90 

0. 

1.86 

3.09 

i 

0. 

0419 

0. 

0350 

4 

3. 

47 

2. 

90 

.1988 


ROLL 


Ang.Vel. rad/sec 

0. 

0.0182 

0.0450 

0. 

0.0453 

0.0745 

0. 

0.0969 

0.1259 

Ang. Accel. rad/sec2 

0/ 

0.0322 

0.0505 

0. 

0.0640 

0.0655 

0. 

0.0785 

0.0935 

Tang. Vel. ft/sec 

0. 

1.35 

3.72 

0. 

3.75 

6.17 

0. 

8.03 

10.42 

Tang. Accel. ft /sec2 

0. 

2.45 

4.18 

0. 

5.30 

5.43 

0. 

6.51 

7.75 

Norm. Accel. ft /sec 2 

0. 

0.02 

0.27 

0. 

0.48 

0.90 

0. 

1.96 

3.09 
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(Continued) 


TABLE  III 
(Continued) 

RESULTS  OF  SHORT  TFRM  POINT  MOTION  PREDICTIONS1 
(RMS  of  maximum  values) 


STATE  4 STATE  5 STATE  6 


Heading  Angle*1 

Slg.Ht.  • 6.5ft. 

180°  135°  90° 

Slg.Ht 

180° 

- io.: 

135° 

> ft. 
90° 

Slg.Ht . - 16, 
180°  135° 

.9  ft. 
90° 

POINT 

4,  SPEED  0 kts. 

VERTICAL  MOTION 

Dbl. Amplitude,  ft. 

2.06 

2.23 

1.52 

3.58 

3.54 

2.52 

5.78 

5.46 

4.24 

Velocity,  ft/sec 

1.96 

2.19 

1.58 

2.98 

3.11 

2.21 

4.14 

4.11 

3.03 

Accel.,  ft/sec* 

2.05 

2.27 

1.85 

2.72 

2.96 

2.25 

3.33 

3.53 

2.59 

TRANSVERSE  ACCELERATION 

Trans .Accel . f t/sec2 

0. 

1.16 

1.72 

0. 

1.49 

2.20 

0. 

1.71 

2.52 

POINT 

4,  SPEED  3 kts. 

VERTICAL  MOTION 

Dbl. Amplitude,  ft. 

2.04 

2.19 

1.56 

3.63 

3.55 

2.55 

5.88 

5.50 

4.26 

Velocity,  ft/sec 

2.12 

2.35 

1.58 

3.35 

3.38 

2.22 

4.68 

4.47 

3.06 

Accel,  ft/sec^ 

2.38 

2.66 

1.80 

3.35 

3.50 

2.23 

4.20 

4.20 

2.61 

TRANSVERSE  ACCELERATION 

Trans. Accel. ft /sec^ 

0. 

1.26 

1.68 

0. 

1.66 

2.18 

0. 

1.95 

2.52 

POINT 

4,  SPEED  5 kts. 

VERTICAL  MOTION 

Dbl .Amplitude,  ft. 

2.00 

2.15 

1.59 

3.63 

3.53 

2.57 

5.91 

5.52 

4.28 

Velocity,  ft/sjc 

2.21 

2.43 

1.60 

3.56 

3.52 

2.25 

5.00 

4.69 

3.08 

Accel.,  ft/sec 

2.61 

2.89 

1.78 

3.78 

3.85 

2.23 

4.78 

4.63 

2.63 

TRANSVERSE  ACCELERATION 

Trans .Accel .ft/s ec^ 

0.0 

1.30 

1.65 

0.0 

1.76 

2.17 

0.0 

2.09 

2.53 

^Using  experimentally  generated  roll  RAO's. 

^Angle  between  direction  of  ship  and  direction  of  waves. 

180°  ■ head  seas,  135°  - bow  seas,  90°  ■ beam  seas. 

•*To  obtain  values  for  Points  1 and  2 multiply  Point  3 values  cited  by  the 
ratio  of  their  vertical  distances  to  the  center  of  gravity,  that  is,  multiply 
by  .840  and  by  .920  respectively. 
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4.0  Design  Criteria  Development 


The  objective  of  this  section  is  to  develop  a set  of  curves 
to  be  used  by  the  ship  designer  in  establishing  the  limits  of 
linear  and  angular  motions  (displacements,  velocities  and  acceler- 
ations) based  on  desired  probabilities  of  exceedance  or  frequencies 
of  occurrence  over  a period  of  ship  operating  time. 

4.1  Current  State-of-the-art 


For  several  years  the  process  of  predicting  wave-induced  loads 
for  design  purposes  has  consisted  of  the  following  steps  and  assump- 
tions. 

1.  It  is  usually  assumed  that  ocean  waves  can  be  represented  by 
"Pierson-Moskowitz"  spectra,  and  that  this  spectral  form  could 
be  considered  sufficiently  "narrow"  to  allow  a Rayleigh  distri- 
bution to  be  used  to  predict  the  amplitudes  of  the  peaks  and 
troughs.  Several  workers  (Reference  3,  for  example)  have 
recently  expressed  doubts  as  to  the  validity  of  these  assump- 
tions, particularly  with  reference  to  the  extreme  conditions 
that  must  generate  design  loads.  The  Pierson-Moskowitz  spectra 
were  Intended  to  represent  fully  developed  sea  waves  (Reference  4). 
Fully-developed  seas  are  only  generated  if  the  wind  blows  for  a 
long  enough  duration  over  an  extensive  enough  fetch.  The 
requirements  for  both  duration  and  fetch  become  larger  as  the 
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wind  speed  increases  and,  for  severe  storms,  become  quite 
unrealistic  as  shown  in  Table  IV  and  Figure  11,  so  anything 
approaching  a "fully-developed"  extreme  sea  state  Is  extremely 
Improbable.  For  example.  Figure  11  shows  that  above  Sea  State 
7 the  storm  duration  required  for  a fully  arisen  sea  exceeds 
the  total  time  expected  for  that  sea  state  in  one  year.  The 
theoretically  attainable  wave  heights  also  shown  in  Figure  11 
will  almost  certainly  not  be  reached  in  the  higher  sea  states 
for  the  same  reason,  so  that  the  highest  waves  and  consequently, 
the  highest  wave-induced  loads  may  actually  occur  in  lower  sea 
states  due  to  their  greater  duration.  In  a similar  way  the 
Rayleigh  distribution  function  is  very  adequate  to  describe 
light  to  moderate  sea  states  but  does  not  describe  severe  sea 
states  with  much  realism.  The  fact  that  the  highest  loads 
are  to  be  expected  in  the  more  frequently  occurring  medium 
sea  states  rather  than  in  the  rare  extreme  sea  states  was 
pointed  out  in  Reference  6,  based  on  full-scale  experimental 
evidence.  There  has  recently  been  a tendency  at  NSRDC  to  pro- 
pose that  calculations  of  extreme  loads  should  be  made  from 
measured  extreme  wave  spectra  but  this  method  is  inconsistent 
with  the  evidence  referred  to  above. 

2.  The  second  major  assumption  is  that  of  linear  response.  Clearly 
if  responses  such  as  heave,  pitch,  bending-moment,  vertical  ac- 
celeration, etc.  can  be  considered  to  be  linear  with  respect  to 
the  wave  then  the  process  of  calculating  extreme  motions  becomes 
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Table  IV  Standard  Sea  State  Parameters  (Reference  5) 


i 


rather  simple.  This  is  a very  real  consideration  and  many 
authorities,  notably  E.  V.  Lewis,  for  example,  in  Reference  7, 
maintain  that  it  is  preferable  to  use  a manageable  design  tech- 
nique that  is  known  to  be  somewhat  approximate  rather  than  a 
more  exact  analysis  that  is  too  complex  to  apply  in  practice 
at  this  time.  Largely  for  this  reason  the  linear  response 
theory  has  been  adhered  to  long  after  it  has  clearly  and 
repeatedly  been  proven  not  to  apply  to  high  performance  marine 
vehicles.  One  reason  for  this  has  been  that,  for  conventional 
ships.  References  6,7,  8,  9 and  many  others  have  shown 

that  the  linear  assumption  gives  very  acceptable  agreement  with 
full-scale  and  model  tank  measurements.  For  planing  craft, 
surface  effect  ships,  etc.,  the  craft  geometry  and  the  dynamic 
forces  acting  on  the  hull  vary  in  a markedly  nonlinear  fashion. 
For  the  case  of  roll  motion  which  is  a highly  nonlinear  mode 
and,  more  specifically  for  rolling  motions  of  a hard  chine  hull 
such  as  that  of  the  T-AGOS,  the  assumption  of  linearity  no 
longer  holds  and  a different  course  of  action  must  be  followed 
in  to  order  to  predict  the  long  term  extremes  for  design  pur- 
poses. For  example,  Figure  12  shows  typical  peak  vertical 
acceleration  distribution  for  the  SES-100B  testcraft  in  Sea 
State  3.  The  distribution  used  here  to  approximate  the  data  is 
the  Weibull  distribution  (to  be  discussed  later)  of  which  the 
Rayleigh  and  Exponential  distributions  are  special  cases  (all 
that  changes  is  the  gradient  of  the  curve) . In  the  case  of 
Figure  12  the  Weibull  distribution  which  fits  the  data  has  a 
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gradient  which  falls  between  that  of  the  Rayleigh  and  Exponential 
as  is  clearly  demonstrated.  Thus  for  this  particular  case  one 
would  use  a distribution  with  the  best  fitting  gradient  to  pre- 
dict the  long-term  extreme  motions  and  loads  given  a certain 
acceptable  probability  or  frequency  of  occurrence  established 
apriori. 

On  the  basis  of  these  observations  it  is  hypothesized  that  the 
Weibull  distribution  may  be  said  to  represent  long  term  roll 
data  in  a reasonable  manner,  and  is  a clear  improvement  on 
either  the  exponential  or  the  Rayleigh  distributions  which  are, 
in  turn,  one-parameter  special  cases  of  the  Weibull  distribu- 
tion. The  Weibull  distribution  is  also  acceptable  from  the 
point  of  view  that  it  is  closely  related  to  the  Rayleigh  dis- 
tribution which  is  the  generally  accepted  representation  of 
wave-heights  and  wave-induced  phenomena. 

In  order  to  make  use  of  the  Weibull  distribution  in  the  analysis 
of  extreme  values  some  of  its  characteristics  muct  be  under- 
stood. These  are  presented  in  Appendix  (I). 

4.2  Prediction  of  long-term  extreme  values 

The  rms  values  predicted  by  the  MIT  5-DOF  program  can  now  be 
used  in  conjunction  with  the  Weibull  distribution  in  order  to  pre- 
dict the  long  term  extreme  values  of  motions,  including  displace- 
ments, velocities  and  accelerations.  Having  obtained  these  long 


term  distribution  curves  then  the  design  criteria  can  be  applied 
in  order  to  obtain  the  desired  limit  motions  and  their  frequency 
of  occurrence  from  the  curves. 

As  discussed  in  Appendix  (I)  the  Weibull  distribution  prob- 
ability density  function  is  given  by: 

f(x)  » c.\C  xc  1 exp  {-(Ax)C}  (1) 

The  value  of  \ can  be  defined  in  terms  of  the  r.m.s.  value  r 
and  the  Gamma  function  T: 


The  cumulative  probability  can  be  expressed  as: 

P(Xj)  - exp  {-(\x)C}  (3) 

or  in  terms  of  r: 

P(x^)  - exp  {-(-^-)C}  (4) 

P(x^)  is  the  probability  of  exceeding  a value  x^  which  in  our 
specific  case  are  amplitudes  of  displacement,  velocity  and  accelera- 
tion. The  exponent  c is  a constant  which  represents  the  gradient 
of  the  distribution  curve.  Clearly  equation  (4)  tells  us  that  de- 
pending on  the  value  of  :he  exponent  c,  the  degree  of  conservativism 
in  the  prediction  of  the  long-term  extremes  will  vary. 

The  effect  of  the  choice  of  distribution  on  the  long  term 
motions  may  be  seen  most  graphically  in  Figure  1.4  of  Appendix  (I). 
At  the  higher  levels  of  probability,  which  are  those  of  interest 
in  establishing  extreme  values,  the  long-term  cumulative 
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distribution  curves  fan  out  over  a wide  range  for  different  values 
of  the  parameter  c.  This  difference  may  be  reduced  somewhat  if  it 
is  hypothesized  that  these  distributions  should  not  be  compared 
with  respect  to  the  rather  abstract  parameter  x£  but  should  rather 
be  compared  on  the  basis  of  the  ratio  x/r  of  x to  the  rms  value. 

The  variation  of  r with  variation  of  the  parameter  c is  shown  in 
Figures  1.1  and  1.2  of  Appendix  ( I ).  Figure  13,  therefore,  shows 
the  long  term  distributions  of  x/r  for  various  c values. 

Now  each  run  in  a test  tank,  of  the  full-scale  vessel  or  a 
computer  run  with  the  motions  program  is  necessarily  a short-term 
event.  During  each  short-term  event,  it  is  assumed  that  each 
ship's  speed  and  heading  and  the  sea  state  remain  constant.  The 
operational  life  of  a ship  can  be  considered  as  a summation  of  a 
very  large  number  of  short  term  events.  If  the  behavior  in  each 
short-term  event  is  known,  and  also  if  the  manner  of  distribution 
of  short-term  events  throughout  the  ship's  life  is  known,  then  a 
long-term  picture  of  the  ship's  life  can  be  built  up.  This  is  the 
method  used  here  to  predict  the  long-term  behavior  of  the  vehicle 
so  that  eventually,  a single  long-term  distribution  of  motions  (or 
loads  or  stresses)  can  be  built  up  to  represent  the  ship's  total 
experience  in  all  speeds,  sea  states,  headings  and  loading  condi- 
tions. In  order  to  accomplish  this,  use  is  made  of  the  description 
of  the  operational  environment  as  shown  in  figure  (11)  as  well  as 
the  operating  envelope  (forward  speed  vs  sea  state) . 


ai 

Weibull  Cumulative  Distribution  of  x/*ma>  (ref.  5) 
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If  the  probability  of  exceeding  a given  motion  x at  a speed 


V,  in  sea  state  S,  at  heading  H and  displacement  W,  is 

then  the  total  probability  P(x)  of  exceeding  the  motion  x is  given 

by  summing  the  component  probabilities  for  all  conditions 


HX)  ' V B I W W VV W>  <S> 

where  P , P , P , P are  the  probabilities  of  occurrence  of  each 
V H 5 W 

velocity,  heading,  sea  state  and  displacement  respectively. 

For  the  purposes  of  this  effort  however,  a few  simplifications 
will  have  to  be  made  since  no  information  has  been  provided  on  the 
operating  envelope  and  the  displacement  probabilities  of  the  T-AGOS. 
Thus  the  combined  probability  of  exceeding  a motion  level  x is 
taken  here  as 


P(x)  - E P_*Pc(x)  (6) 

S S 5 

where  Py,  P^,  and  P^  are  not  taken  into  consideration.  Instead, 
P(x)  is  computed  using  equation  (6  ) for  specified  displace- 
ment, heading  and  speed.  However,  as  the  information  on  the 
operating  envelope  is  made  available,  then  a more  rigorous  calcu- 
lation could  be  conducted.  The  use  of  equation  (6)  will  yield  a 
conservative  estimate  since  P^,  , and  Pw  are  not  considered. 

5.0  Design  Criteria  Recommendations 

A reasonable  criteria  for  selecting  the  design  motion  level 
for  the  T-AGOS  is  to  select  that  displacement,  velocity  or 
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acceleration  which  would  occur  once  In  a certain  period  of  opera- 
tion. For  example,  based  on  a nominal  2000  hours  per  year  at 
sea,  the  frequency  of  occurrence  F,  of  once  In  2000  hours  can  be 
expressed  as 

T ■ 2000  hours  x 3600  seconds 

hour 

■ 7.2  x 10^  secs 

F - 1/T  - 1 , - 0.139  x 10-6 (sec-1) 

7.2xl06 

Now  for  any  single  short-term  distribution  of  the  frequency  of 
occurrence  F (x)  of  any  parameter  x is  given  by 

F (x)  - nP (x)  (7) 

where  n is  the  number  of  occurrences  per  unit  time  and  P(x)  is  the 
probability  of  exceeding  x in  any  one  excursion.  Both  f(x)  and  n 
have  dimensions  of,  for  example,  (sec  1) . Thus  the  combined  fre- 
quency F (x)  of  exceeding  any  particular  level,  x. of  motion  response 
is  obtained  from  equations  (6)  and  (7)  and  is  given  by 

F (x)  - nP(x)  - Z n Pg-P^x)  (8) 

S 

The  parameter  n is  a function  of  displacement,  speed,  sea  state 
and  heading. 
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The  average  period  of  rolling  is 


I (9) 

where  ri^  ■ rms  roll  amplitude 

n4  ■ rms  roll  velocity 

The  number  of  occurrences,  n,  per  unit  time  is  then  given  by 

n " T (10) 

Tn  2TTT14 

u 

where  04  and  n4  are  the  rms  values  for  given  sea  state,  heading, 
displacement  and  speed. 

The  value  of  n given  by  0-0)  is  then  substituted  in  (8)  in 

order  to  obtain  the  combined  frequency  of  occurrence  F(x). 

For  our  purposes  we  will  consider  the  beam  seas  roll  case  as 
being  the  most  critical.  This  is  indicated  by  the  output  of  the 
computer  simulations  where  the  largest  rms  amplitudes  of  roll 
occurred  in  beam  seas  and  increased  with  increasing  sea  state. 

For  the  beam  sea  condition  forward  speed  (0,  3,  6 knots)  had  no 
effect  on  roll. 

Table  (V)  summarizes  the  data  needed  to  develop  long  term 

distribution  curves  for  each  of  the  six  parameters  of  interest, 

i.e.,  p. , n, , p. » V . a , and  a . In  addition  to  the  rms  values 
4 4 4 c t n 

of  these  variables  as  a function  of  sea  state  for  beam  seas  and 
forward  speeds  of  0,  3,  and  6 knots,  Table  (V  ) also  shows  the 


probability  of  occurrence  of  each  sea  state,  P^,  which  has  been 
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TABLE  V Roll  motions  as  a function  of  Sea  State  Of  = 90  ; V = 0,  3,  5 knots;  R = 82.83') 


calculated  from  figure  11.  Likewise  the  value  of  n and  T are 
given  for  each  sea  state.  In  all  cases  the  value  of  R (the  dis- 
tance from  the  c.g  to  the  point  In  question)  Is  taken  to  be  the 
maximum  of  82.83  ft. 

Since  at  the  present  time  there  Is  no  experimental  data 

(model  or  full-scale)  which  would  allow  one  to  determine  the 

value  of  the  constant  c in  the  Weibull  distribution,  it  was  felt 

that  a reasonable  path  to  follow  here  was  to  present  long  term 

distribution  curves  for  the  six  parameters  of  interest  using  values 

of  c of  1.3  and  2.0.  With  c ■ 1.3  the  design  levels  would  be 

more  conservative  than  in  the  case  of  c ■ 2.0  which  corresponds 

to  the  Rayleigh  distribution.  This  is  clearly  illustrated  in 

figure  (13).  For  example,  the  probability  of  exceeding  the  ratio 

x/x  ■ 3.0  is  .00001  when  c ■ 2.0  while  the  same  ratio  is 
rms 

exceeded  with  a probability  of  .006  when  c ■ 1.3.  In  this  effort 
the  value  of  c ■ 1.3  has  been  chosen  based  on  observations  and 
analyses  of  other  non-linear  phenomena  (refs.  5,  10,  11.  12). 

Figures  (14)  through  (19)  show  the  results  of  this  procedure 
for  all  six  variables.  In  each  case  design  values  have  been 
obtained  applying  the  criteria  of  one  exceedance  per  2000  hours. 
Conservative  and  less  conservative  estimates  are  presented  depend- 
ing on  the  choice  of  distribution.  The  former  corresponds  to  a 
Weibull  distribution  with  parameter  c - 1.3  while  the  latter  is 
a Weibull  distribution  with  c ■ 2.0  (i.e.  the  Rayleigh  distribu- 
tion). It  is  pointed  out  that  these  estimates  already  have  a 
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built-in  conservativism  due  to  the  assumption  of  2000  hours  of 
operation  in  beam  seas  which  is  not  operationally  realistic.  To 
avoid  this  situation  one  must  include  the  probability  of  heading 
occurrence  in  the  computation  of  the  total  combined  probability 
or  frequency  of  exceedance  according  to  equation  (5).  However, 
the  curves  of  figures  (14)  through  (19)  do  take  into  account  the 
expected  frequency  of  occurrence  of  all  sea  states. 

In  view  of  the  above  assumptions,  and  the  limitations  caused 
by  not  having  an  operational  envelope  defined  prior  to  conducting 
this  study,  the  long-term  design  values  presented  here  should  be 
used  with  caution  until  a more  complete  analysis  is  conducted.  As 
indicated  earlier  there  is  no  long-term  experimental  roll  data  for 
the  T-AGOS  or  similar  hulls.  Until  such  information  becomes  avail- 
able a value  of  c ■ 1.3  would  probably  be  satisfactory  and  conserva- 
tive compared  to  a choice  of  c * 2.0.  However,  one  must  consider 
the  added  degree  of  built-in  conservativism  associated  with  the 
heading  assumption  (2000  hours  in  beam  seas).  Consequently,  in 
order  to  compensate  for  this  factor,  it  is  recommended  that  for  now 
the  less  conservative  design  value  be  used.  It  is  pointed  out  that 
the  use  of  the  criteria  "one  exceedance  in  2000  hours"  is  simply  an 
example  to  illustrate  the  use  of  the  long  term  distribution  curves. 
The  choice  of  the  actual  criteria  to  be  applied  is  the  responsi- 
bility of  the  SHAPM. 

Finally,  it  would  appear  to  be  very  desirable  to  devote  some 
effort  to  estimating  with  some  degree  of  confidence  the  value  of 
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the  parameter  c.  A possible  way  of  doing  this  is  to  analyze  the 
tapes  obtained  from  past  and/or  future  tow  tank  tests  where  roll 
motion  has  been  measured.  With  this  information  available  the 
peaks  and  troughs  of  the  responses  could  be  counted,  and  short- 
and  long-term  distribution  plots  similar  to  those  of  figure  (12) 
could  be  obtained.  This  would  allow  one  to  determine  the  value  of 
c which  best  fits  the  experimental  data  and,  consequently,  the  one 
to  use  in  developing  long-term  distributions  for  use  in  design. 
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Appendix  I 


CHARACTERISTICS  OF  THE  WEIBULL  DISTRIBUTION 


The  generalized  gamma  distribution  is  described  by  the  probability 
density  function 


f(x) 


c .cm  cm-1  / ,,  .c-i 
A x exp  i-(Xx)  } 


T(m) 


where  c,  X and  m are  the  three  parameters  of  the  distribution. 


The  cumulative  density  function  is  given  by 
x 

f (x)  dx 


Q(x) 


•/ 


note  asx  + « 


Q(x) 


/ 


f (x)  dx 


r(m) 


cm 


cm-1 


exp[-(Xx)cJ  dx 


c 

r<m) 


.cm 


£iki 

ck 


where  k ■ m in  this  case 


cX 


- 1 
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(1) 


(2) 
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The  mean  value  of  the  distribution  is  given  by 
Tt  ■ ^f(x)  * x dx 

flO 

“ xcm  exp[-(Ax)d]  dx 


_c  cm  rot l 
^ cA* 


cA”1  T(m  + ^ ) 


where  k ■ m + — 
c 


xA  - 


r(m)  cA 

r(m  + i ) 

r(o)A 


T(m  + “ ) 


cm+1 


and 


r(m) 


(3) 


(3a) 


The  standard  deviation  o about  the  mean  value  is  given  by 

O - [ /f(x)(x  - x)2  d x|1/2  ( 

' o 

• j'ffe'  xCm\j  xCmfl  exP  t-(**)C]  dx 

o 

OS  os 

J xcm  exp  [-(Ax) c]  dx  + TfZj  xCWrml  exp[-(Ax)c]  dxjj  1/2 


o 

C .Ctl 


T(m) 


- j£l2H_± 

I r(m) AJ 


2/cl  > -=2  U/2 


and 


oA  ■ . (yX)2ll/2 


r(m) 


I 


48(a) 


The  root  mean  square  value  abouC  Che  x * o axis  Is  given  by 


•1/ 


f(x)  x*  dx 

1/2 


1/2 


(5) 


TCm  + 2/c) 

r(m)A2  j 

r(m  + 2/c)  j1^2 
r(m)  j 

2 — 2 1/2 

Note  that  • (o  + x ) 
mis 


and 


W “ 


The  values  of  Ax,  Ao  and  Ax^„  are  plotted  for  a range  of  values  of  c and 

rms 

m In  Figure  1.1. 


It  is  pointed  out  that  the  generalized  gamma  distribution  reduces 
to  the  Welbull  distribution  in  the  special  case  when  m ■ 1 and 
that  the  Weibull  distribution  reduces  to  the  exponential  distri- 
bution when  c ■ 1 and  to  the  Rayleigh  distribution  when  c - 2. 

These  three  distributions  are  indicated  on  Figure  1.1. 

The  Welbull  probability  density  function  is  obtained  by  substituting 
m ■ 1 in  equation  1 which  leaves  A and  c as  the  two  parameters: 

f(x)  ■ cAc  xc_1  expl-(Ax)c]  (6) 

The  mean,  mode,  median  and  rms  values  are  given  by 

Ale  - T(1  + 1/c) 

’ [<c-l)/c]l/c 

■ I-Ml))1'' 

and 

‘-rma  ’ [r(l«/c)]l/2 

from  equations  3(a)  and  5(a). 
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The  cumulative  probability  function  is  given  by 


P(x)  - 1 - Q(x) 


■f 


f(x)  dx 


exp  (-(Ax)c] 


(7) 


f(x)  and  P(x)  are  plotted  for  a range  of  values  of  c in  Figures  1.2, 
1.3  and  1.4.  Other  terms  that  are  often  used  to  describe  statisti- 
cal distributions  are  the  "average  of  the  highest  one-third,  x^" 
and  the  "average  of  the  highest  one-tenth,  x1/10"  etc.  In  general, 
the  "average  of  the  x^"  values  as  defined  in  Figure  1.5  can  be 
evaluated  as  follows. 


The  value,  x^^,  a^ove  which  1/b  of  the  values  of  x are  expected 
to  fall  is  given  by  putting  P(x)  - 1/b  in  equation  7 and  solving 
Cot  x1/b: 


i.e. 


b “ exp  ^~^xi/b^  ^ 
Xx1/b  - {-in(l/b)}1/c 


(8) 


The  average  Ax^yb  of  the  highest  1/b  values  is  therefore  given  by 


f 


x . f(x)  dx 


At,  „ - *l/b 


1/b 


f 

^1/b 

{ 

1/b 


b I cAc  xC  exp  l-(Ax)C]  dx 


(9) 
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This  can  be  evaluated  numerically.  The  probability  Z(x^^)  of 


exceeding  the  value  x^b  is  given  by: 


p(*l/b>  ■ exp  {-(Xx1/,b)C) 


which  can  also  be  evaluated  numerically. 


Values  of  x^b  and  are  plotted  in  Figures  1.6  and  1.7 


respectively.  In  Figure  1.7,  it  may  be  noted  that  the  average 
of  the  highest  1/3  and  1/10  values  do  not  vary  greatly  from  one 
value  of  c to  another. 
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Characteristics  of  *-he  Weibull  Distribution. 


